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: ABSTRACT 

We present a detailed investigation of the young stellar populations (YSP) in 
O ■ the radio-loud ultra luminous infrared galaxy (ULIRG) PKS1345+12(z=0.12), 



based on high resolution Hubble Space Telescope (HST) imaging and long 
slit spectra taken with the William Herschel Telescope (WHT) on La Palma. 
While the images clearly show bright knots suggestive of super star clus- 
ters(SSC), the spectra reveal the presence of YSP in the diffuse light across the 



full extent of the halo of the merging-double nucleus system. Spectral synthe- 
sis modelling has been used to estimate the ages of the YSP for both the SSC 
and the diffuse light sampled by the spectra. For the SSC we find ages tssc<6 
Myr with reddenings 0.2 < E(B - V) < 0.5 and masses 10 6 < M^§g < 10 7 
Mq. In the region to the south of the western nucleus that contains the SSC 
our modelling of the spectrum of the diffuse light is also consistent with a 
relatively young age for the YSP (~5 Myr), although older YSP ages can- 
not be ruled out. However, in other regions of the galaxy we find that the 
spectra of the diffuse light component can only be modelled with a relatively 
old post-starburst YSP (0.04 - 1.0 Gyr) or with a disk galaxy template spec- 
trum. The results demonstrate the importance of accounting for reddening 
in photometric studies of SSC, and highlight the dangers of focussing on the 
highest surface brightness regions when trying to obtain a general impression 
of the star formation activity in the host galaxies of ULIRGs. The case of 
PKS 1345+1 2 provides clear evidence that the star formation histories of the 
YSP in ULIRGs are complex. While the SSC represent the vigorous phase 
of star formation associated with the final stages of the merger, the YSP in 
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the diffuse fight are likely to represent star formation in one or more of the 
merging galaxies at an earlier stage or prior to the start of the merger. 

Intriguingly, our long-slit spectra show line splitting at the locations of 
the SSC, indicating that they are moving at up to 450 km s _1 with respect to 
the local ambient gas. Given their kinematics, it is plausible that the SSC have 
been formed cither in fast moving gas streams / tidal tails that are falling back 
into the nuclear regions as part of the merger process, or as a consequence 
of jet-induced star formation linked to the extended, diffuse radio emission 
detected in the halo of the galaxy. 

Key words: Radio galaxies, ULlRGs, Young stellar populaitons, mergers, 
starbursts. 



1 INTRODUCTION 



A variety of surveys over the last two decades have re vealed populat i ons of 



the bulk of their r a diation at infrared waveleng t hs ( Soifer fc et a] 



198 



4 



Houck et al. 



1985 



Le Floc'h et al. 



1984a 



2005; iPerez-Gonzalez et al. 



ga. 



b 



axies that emit 



Houck &: et all . 



20051 ). Such objects 



are thought to be resp onsible for most of the star formation activity in the distant Universe 
( iLe Floc'h et all 120051 ). while in the local Universe these galaxies are classified as Luminous 
(L ir > 10 11 L ) or Ultraluminous (L ir > 10 12 L ) infrared galaxies (LIRGs/ULIRGs). The 
morphologies of ULIRGs clearly suggest that most of them are recent or ongoing mergers. 
The presence of an AGN in a subset of these objects is also well known, with some ULIRGs 
classified as QSOs a t optical wavelengths. For a more complete review of the properties of 
LIRGs/ULIRGs see lsanders fc Mirabell ( 1996 ). 

Much attention has been paid recently to the so-called warm ULIRGs (mid-infrared 
colours of /25//60> 0.2)0. These objects represent the ~20-25% of the total ULIRG popula- 
tion discovered by the IRAS satellite. Most of them have AGN optical spectra, very large 
molecular gas masses (Mm,, ~ 10 10 M^: Sander s et al. 1988) and are found in an advanced 
merger state. In fact, ISurace fc Sanders! (Il999l ) found that 75% of the warm ULIRGs in 
their sample of 12 objects were associated with single nucleus galaxies. These properties 
suggest that suc h objects represent a transitional stage between cool ULIRGs and Radio 
Galaxies/QSOs ( ISanders et all Il988l ): cool ULIRGs evolve into warm ULIRGs on their way 
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1 The quantities /2s and /go represent the IRAS flux densities (non-colour corrected) in units of Jy at 25 and 60 fim respectively 
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to becoming QSOs/radio galaxies. Indeed, models predict that the tidal forces associated 
wit h mergers lead to large concentrations of gas and dust in the nuclear regions of the galax- 



ies (jMihos fc Hernquist 



1996 



Barnes fc Hernquistl . Il996l ); such concentrations may trigger 



both AGN and starburst activity. However, the timing of the AGN relative to the major 
starburst triggered by the merger remains uncertain. Although some models assume that 
both starburst and AGN will occur simultaneously, w hen there is sufficient concentration 
of gas in the nuclear regions (Idi Matteo et all 120051 ). the idea that cool ULIRGs evolve 
eventually into QSOs (jSanders et all Il988l ) suggests that there may be a significant delay 
between the starburst and visible AGN activity. 

Given that the models make specific predictions about the histories of the star formation 
triggered in the course of major gas-rich mergers, studies of the stellar populations in ULIRGs 
provide useful information about the mergers and, potentially, allow us to test the models. 
Moreover, for those ULIRGs with AGN and powerful radio jets, the ages of the young 
stellar populations (YSP) can be used to establish the order of events during the merger, for 



example, whether t 



re AGN was activated before, at the same t ime as, or after the merger- 



induced starburst (jTadhunter et al. 



2005 



Emontsetal 



20061 ). Thus studies of the YSP 



have the potential to help us to understand the possible links between ULIRGs and AGN. 

In this context, there have been surprisingly few studies of the YSP in ULIRGs. Two 
complementary approaches have been used in the past: one based on photometric analysis of 



images of the galaxies taken at different wavelengths, 



so dated with the bright kno ts detected in the obje cts (jSurace et al 



1999 



Surace fc et all . 



2000£ 



bcused mainly on studing the YSP as- 



1998 



Wilson et al. 



Surace fc Sanders! . 



20061 ) ; the other comprising spectrocopic anal y 



SIS o 



2001 



2000 



b, 



JV/optical spectra of m ore spatially extended regions ( ICanalizo h Stockton 
Tadhunter et all I2OO5L Note that such optical studies sample the star formation in 



the systems prior to the current prodigious formation activity represented by the ULIRG 
that is likely to be hidden from view at optical wavelengths. In this paper we combine op- 
tical spectroscopy and imaging observation in order to investigate the YSP in the warm 
ULIRG /Radio Galaxy PKS1345+12 



(Evans et al.. 


1999: 


Holt et al.. 


2003) 



in an elliptical-like host galaxy (e.g. Axon et al. 2000). The galaxy has a clear double nucleus 
separated by ~1.8 arcsec or 4.3 kpfl and the AGN activity is associated with the western 



2 Ho = 75 km s 1 , qo = 0.0 assumed throughout resulting in a scale of 2.37 kpc arcsec 1 at z = 0.122 
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nucleus, which has a spectrum characteristic of narrow line radio galaxies (jEvans et al. 



1999 



Holt et al. 



2003|). 



Si gnificant YSP have also been detected in PKS1345+12 at optical wavelengths (ISurace et al. 



1998 



Tadhunter et al. 



20051 ). which may have been formed as a consequence of the merger. 



Moreover, evidence for prodigious hidden star formation activity is provided by the detec- 
tion of a large far-IR excess by the IRAS satellite (Lm = 1.7 x 10 12 M , (Evans et al. 
1999)). In addition, the CO observations indicate the presence of a substantial reservoir 
(3.3 x 10 10 M ) of molecular gas within ~2.5 kpc of the active (western) nucleus (p > 
2000 M Q pc~ 3 : E yans et al. 1999), c onsist ent with measurements in other ultraluminous 
infrared galaxies (IBryant fc Scovilld . I1999D . The double nucleus, the distorted large-scale 
morphology, the presence of a rich inter-stellar medium (ISM) and YSPs clearly indicate 



that PK S1345+12 represents the later stages o 



galaxy ( IHeckman et al 



1986 



Surace et al 



1998|) 



a merger involving at least one gas-rich 



For the imaging part of the project we have used HST archive data taken with various 
cameras and filters sensitive from the UV to the near-IR. The wide spectral coverage of 
the observations allows us to make accurate estimates of the ages of the stellar population 
associated with the bright knots identified in the images. Fo r the spectroscopy we have used 



high-quality long-slit spectra presented in 



Holt et al. 



(120031 ) . The spectroscopic data enable 



us to investigate the ages and mass contributions of the YSP detected at optical wavelengths 
in the diffuse halo of the galaxy. We compare the results with model predictions in order to 
understand the past and future of PKS1345+12. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 HST DATA 

The HST dataset comprises four sets of images taken with the FOC, WFPC2, ACS and 
NICMOS cameras. A summary of the observations is presented in Table 1. 
The observations and data reduction of the FOC images are described in 



Hurt et al. 



(119991 ). These pre-costar images were taken with the f/96 relay, with a combination of the 
F320W filter and three different polarizing filters (POLO, POL60 and POL120). The central 
wavelength of this configuration is A c = 3100 A for a power-law spectrum /\ocA -1 . 

For the purpose of this paper we are not interested in the polarization, but rather the 
total flux at various wavelengths. Therefore we adapted the data to make them suitable for 
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Camera 


Filter 


A c (A) (rest frame) 


AA(A)(rest frame) 


Exposure Time (sec) 


FOC 


F320W-POL0 


2766 


752 


602 


FOC 


F320W-POL60 


2773 


752 


529 


FOC 


F320W-POL120 


2766 


752 


602 


WFPC2 


F439W 


3844 


423 


2200 


WFPC2 


F814W 


7127 


1567 


835 


ACS(HRC) 


FR459M 


4539 


410 


2480 


ACS(HRC) 


F550M 


4973 


488 


2480 


ACS(WFC) 


FR647M 


5897 


530 


800 


ACS(WFC) 


FR647M 


6549 


588 


800 


NICMOS 


F110W 


10060 


5349 


95 


NICMOS 


F160W 


14316 


3566 


95 



Table 1. Summary of the HST observations of PKS1345+12 used in this paper 



knot 








flux(/10" 


-18 er g gec 1 


cm" 2 A -1 ) 








number 


F320W 


F439W 


FR459M 


F550M 6 


FR647M 


F647M 6 


F814W 


F110W 


F160W 


CI 


2.86±0.57 


2.09±0.15 


1.55±0.06 


1.35±0.05 


0.75±0.05 


1.30±0.05 


0.70±0.05 


<0.40 


<0.24 


C2 


1.44±0.29 


1.68±0.12 


1.35±0.05 


1.11±0.04 


0.98±0.05 


1.43±0.06 


0.69±0.05 


<0.40 


<0.24 


C3 a 










0.88±0.05 


1.35±0.05 


0.45±0.04 






C4 


1.56±0.30 


1.17±0.10 


0.71±0.03 


0.76±0.03 


0.58±0.04 


0.96±0.04 


0.35±0.04 


<0.40 


<0.20 



Table 2. Fluxes measured for the clusters at various wavelengths from the HST images. 

"For C3 we were only able to obtain accurate photometry for three photmetric bands, due to the proximity 

of this cluster to one of the nuclei (see Fig 1). 

'These two filters contain the [Om] and Ha emission lines. 

this study. Since the three images were taken at the same wavelength, we averaged them 
using the IRAF routine IMCOMBINE. Assuming that the light from the bright knots is 
unpolarized — a reasonable assumption — the polarizing filters will transmit half of the total 
flux. Therefore the fluxes derived from these data using the standard photometric calibration 
factors were multiplied by a factor of two in order to work out the total fluxes. The absolute 



photometric accuracy of the f/96 relay is 10 - 20 % iBauml (119941 ) 



The WFPC2 observations and their reduction are described in detail by 



Surace et al 



( 119981 ). The PKS1345+12 images were centred in the PC chip (0.046 arcseconds pixel x ). 



The filters used were F439W and F814W , which are very similar to the standard Johnson B 



and Cousins I filters (IHoltzman fc et al. 



1995a 



bj). Dark subtraction, bias subtraction and 



flat fielding were carried out using the standard data reduction pipeline procedures at STScI. 
The images were combined using the IRAF/STSDAS routine CRREJ which removes cosmic 
rays in the combined image. The calculated errors on the photometric points, including both 
the photon noise and calibration uncertainity, are ±7-10 per cent. 

Four ACS images of PKS1345+12 are used in this paper. The Wide Field Channel (WFC: 
0.049 arcseonds pixel -1 ) was used in combination with FR647M medium band ramp filter to 
take two of the images, one centred on the redshifted Ha emission line wavelength (hereafter 
Ha image), and the other centred on the nearby continuum (hereafter Ha continuum image). 
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N 







Figure 1. The ACS Ha continuum image of PKS1345+12 showing the presence of two nuclei, a distorted halo and four bright 
compact clusters embeded in this halo. The image also shows the presence of dust features spread to the south of the galaxy. 
These dust features intercept the positions of three of the clusters (C1,C2,C4). 



The High Resolution Channel (HRC: 0.027 arcseconds pixel -1 ) was used in combination with 
F550M medium band filter (A = 5579 A A\ FWHM = 547 A) and FR459M medium band 
ramp filter to take the other images. In the former, the aim was to take an image of the 
galaxy at the wavelength of the [O ill] A5007 emission line (hereafter O ill image), whereas in 
the other images the pivot wavelength was chosen to take an image of the continuum near 
the [Om]A5007 line ( hereafter Om continuum image). The FR647M(A C = 6615 A) image 
is shown in Figure 1. The data were reduced following the standard data reduction pipeline 
procedures at STScI which employ two packages: CALACS package, which performs dark 
substaction, bias substraction and flat field corrections producing calibrated images, and 
MULTIDRIZZLE package which corrects for distortion and performs cosmic rays rejections. 
Since two of the ACS images were set up to measure the Ha and [Om]A5007 emission 
lines, they were not suitable to be used for modelling the continuum, and only the FR647M 
(A = 6615 A) and FR459M images were used for that purpose. The calculated errors of 
the photometric points including both the photon noise and calibration uncertanty are ±5-7 
per cent. 

For the FOC, WFPC2 and ACS data, the routine PHOT in IRAF was used measure 
sky-subtracted fluxes for the clusters within circular apertures typically 3-10 pixels in radius, 
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correcting for aperture losses using a synthetic PSFs generated with the TINYTIM program. 
The measured fluxes are summarised in Table 2. 

The observations and reduction of the NICMOS data are described in Scoville et al 



((20001) ■ The images were taken using camera 2 of NICMOS, ha ving a spatial resolutio n 
of 0.0762 and 0.0755 arcseconds pixel -1 in x and y respectively ([Thompson et all I1998T ). 
The F110W(1.10 £tm, AX FWHM ~ 0.6 /mi), F160W(1.60/mi, AX FWHM ~ 0.4 /mi) filters 
were used for the observations. The images we used for this paper were fully reduced and 
calibrated using the standard pipeline procedure. 

No SSCs are clearly visible in the NICMOS images. However, the images are useful for 
measuring upper limits, which provide constraints on the long wavelength SEDs. For that 
purpose, the TINYTIM program was used to create synthetic point source images for each 
filter which were scaled using a range of flux scaling factors. The synthetic point sources 
were then added to the images at the expected locations of the SSC, and the flux scale 
factor increased until the point sources were just detected in visual inspection of the images. 
The fluxes of the just-detected point sources were taken as upper limits on the true near-IR 
fluxes of the SSC. We summarize the results based on the NICMOS data in Table 2. 



2.2 Spectra 



In addition, we have used long-slit spectra to study the YSP associated with the diffuse 



lig 
in 



it compon e nt in PKS1345+12. The observations and reduction processes are described 



Holt et al 



(120031 ). The spectra were taken in 2001 May with the ISIS spectrograph on 
the 4.2m William Herschel Telescope (WHT) on La Palma. The final wavelength range 
is 3275-6813 A in the blue and 6212-7720 A in the red. Useful spectra were taken along 
two slit positions: PA 160 and PA 230 with a 1.3 arcsecond slit. The final wavelength 
calibration accuracy was 0.059 A and 0.112 A in the red and blue respectively, the spectral 
resolution is 3.66±0.09 A in the red and 4.54±0.10A in the blue, and the relative flux 
calibration acc uracy is ±5 per cent. A preliminary analysis of this dataset for one aperture 



is presented in 



Tadhunter et al 



(120051 ). whereas for the current study several apertures were 



extracted and analysed from both slit positions using the STARLINK packages FIGARO 



and DIPSO. Prior to mod e 
V) = 0.034: Ischleeel et al. 



ling, 



the spectra were corrected for galactic reddening (E(B — 
(119981 )) and we also subtracted the nebular c ontinu um for those 



apertures with strong emission lines (IDickson et al. 



1995 



Holt et al. 



2003 ). In the case 
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of the nuclear aperture, where we have a good estimate of the redd e ning of the various 



kinematic components, the optimum reddening model of iHolt et al.l (120031 ) was used to 
generate a nebular continuum. However, for apertures 1 and 2 the reddening is not known 
accurately because of the effect of the underlying absorption lines on the high order Balmer 
lines. For those apertures we used Ha to generate the nebular continuum and considered 
two extreme cases: (i) maximum nebular continuum assuming no reddening of the emission 
line region; (ii) zero nebular continuum corresponding to high reddening. In the remaining 
apertures (PA 230, A & B; PA 160, 3 & 4) no nebular continuum subtraction was deemed 
necessary. We then performed spectral synthesis modelling of all the extracted spectra; the 
results are presented in next section. 



3 RESULTS 

3.1 Photometric Analysis 

The deep ACS Ha continuum image presented in Figure 1 shows the presence of the two 
nuclei of this double nucleus system, as well as a common non-elliptical, distorted, halo 
of diffuse light surrounding the galaxy. Four bright super are clearly 

identified in the image and are the subject of our photometric analysis. Dust lanes are 
also detected as fine, extended features of lower surface brightness in the extended halo to 
the south of the galaxy. Significantly, the most prominent of these dust lanes intercept the 
positions of three of the SSC (CI, C2, C4). We further note that one of the clusters (C3) is 
close to the western nucleus. Due to this proximity, the flux measurements are less reliable 
for this cluster because of the large gradient in the galaxy background light. Moreover, for 
some of the images this cluster is not clearly detected. Therefore, it has not been considered 
in the modelling process. 

Rather than the colour-colour or colour-magnitude diagrams that have been used previ- 
ously to determine the ag e s of SSC in ULIRGs from multi- wavelength HST imaging observa- 



tions (ISurace et al. 



1998 



Wilson et al. 



20061 ). we instead adopt the approach of modelling 



the optical/UV spectral energy distributions (SEDs) derived from our photometric mea- 
surements. The main advantage of this approach is that it utilises all of the photometric 
information simultaneously, and thereby allows both the reddening and age of each cluster 



3 Althought these bright knots appear unresolved, at the distance of PKS1345+12 we cannot rule out the idea that they 
represent multiple star cluster systems rather than single SSC. 
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to be determined, without the age/reddening degeneracy that can be a problem for single 
colour-colour and colour-magnitude diagrams. 



Stellar population synthesis models in the GALAXEV (version 2003) library (IBruzual fc Chariot . 



20031 ) were used to estimate the ages and masses of clusters CI, C2 and C4 from their re- 



spective SEDs. We used the galaxevpl program within the library to construct synthetic 
SSC spectra with ages in the range lMyr to llMyr. These spectra were c reated as s uming 
an instantaneous burst of star formation, solar metallicity and using the ISalpeterl (119551 ) 
IMF with lower and upper mass cutoffs = 0.1 M and mjj = 100 M Q . Reddened model 
spectra were created for t he reddening; r ange 0.0 < E(B — V) < 1.0 using the parametrized 



Galactic extinction law of 



Seatonl (119791 ). The five photometric UV/optical points suitable 



for the fit, i.e. those not affected by emission line contamination, were modelled using a 
minimum \ 2 technique to find the ages of the SSC. Table 3 presents ranges of ages and 
reddenings for which acceptable fits were obtained (Xred ~ -Qi the ranges of measured masses 
are also presented. During the fitting process, the same (±10%) error was assumed for all 
five photometric points. Taking in account the calibration error and the photon noise, both 
already discussed in section 2.1, plus small shifts in the central wavelengths of the filters, 
and uncertainties due to the gradient in the background near the clusters, we believe that a 
mean error of ±10% is realistic. 

The best fitting models for the three clusters are shown in Figure 2. For comparison, some 
models that fall outside of the range of acceptable fits are also shown. The latter represent the 
models with age/reddening combinations that give the best reduced chi-squared just outside 
the acceptable zone. Taking the example of cluster C2 (see Figure 2b) the fit obtained for 
the 4 Myr with E(B-V)=0.5 model is noticably better than that obtained for either the 
3 Myr with E(B-V)=0.6 or 6 Myr with E(B-V)=0.4 models. Moreover xled va hies found for 
YSP ages of 4-5 Myr and reddening of E(B-V)=0.5 are <1, while those found outside this 
range of ages are > 2. Therefore, the range of ages considered to be valid for the case of C2 
is 4-5 Myr, with the best fit found for a template of 4 Myr with a reddening of E(B-V)=0.5. 
Note that the high fluxes observed for the FOC F320W point for the three clusters, rule out 
starburst ages above 50 Myr, even taking in account any reddening effect, since the observed 
U — B is already bluer than such a starburst. 

Further evidence for the young ages measured for the SSC in PKS1345±12 is provided by 
the detection of emission lines from the clusters in both the HST images (Figure 2) and the 
spectra (see next section). Based on the fluxes determined from our Ha emission line images, 
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Table 3. The results of the spectral synthesis modelling applied to the optical/UV photometric points of Table 2. 



knot 


Age(Myr) 


E(B-V) 


Mass(10 6 M ) 


number 








CI 


<7 


0.2-0.5 


~2.5 - 9.0 


C2 


1-5 


0.5 


~4.0 - 4.8 


C4 


<7 


0.2-0.5 


~1.3 - 5.0 



and correcting for both the SSC continuum and a 40% contribution from [Nil] AA6548, 6584 
emission lines in the filter bandpass, we determine Ha equivalent widths of 250±37, 370±34 
and 606±86 A for clu sters CI, C2 and C4 respectively. The instantaneous burst models of 



Leitherer et al. 



(119991 ) show that such high equivalent widths are only attained for clusters 
with ages t c i < 6 Myr, regardless of the assumed IMF. Since the ISM in the clusters may not 
absorb all the ionizing photons generated by the OB stars (as assumed by the Leitherer et al. 
1999 models), the ages could be significantly less than this upper limit. Thus the detection of 
Ha emission lines with high equivalent width is consistent with the range of ages determined 
from SED modelling of the clustergj. Finally we note that the NICMOS upper limits rule 
out ages older than ~11 Myr with reddenings similar to or lower than those presented in 
Table 3, or ages <7 Myr with higher reddening values than the models in Table 3. 

The masses of the SSC detected in PKS1345+12 (10 6 < M ssc < 10 7 M ) are comparable 
with those of the most massive clusters detected i n the Milky Way (e .g. Omega Centauri), 



m merger 



system s such as the ULIRG Arp220 flWilson et al. 



20061) and Th e Antennae 



( IFall et all 120051 ) . and in some other interacting systems (Ide Grijs et all 120031 ). Given the 
difficulty of detecting individual SSC at the relatively high redshift of PKS1345+12, it is not 
surprising that the clusters we have detected are relatively young and massive. Certainly, 
our imaging observations do not rule out the presence of a more numerous cluster population 
in the halo of PKS1345+12, comprising SSC that are older and/or less massive than CI, C2 
and C4. 



3.2 Analysis of the Spectra 

Figure 3 shows the slit positions and the apertures used for the spectroscopic analysis. The 
same routine as described above was used to create simple stellar population templates 
of different age and reddening. Reddened model spectra were created using E(B-V) values 

4 Note that, although we have not bee able to make detailed fits to the continuum SED of cluster C3, it is clear from Table 1 
that this cluster also has Ho detected at high equivalent width, and therefore it is also relatively young (t c ; < 6 Myr). 
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Figure 2a 




4000 5000 6000 

Rest wavelength (A) 



Figure 2b 




4000 5000 

Rest wavelength (A) 



Figure 2c 




4000 5000 

Rest wavelength (A) 



Figure 2. The results of fitting different templates to the fluxes measured for the three clusters used in the photometric 
analysis (flux measurements shown by solid symbols). Figure 2a shows the fits obtained for Cf using templates of 2 Myr with 
E(B — V) = 0.5 (solid line) and 8 Myr with E(B — V) = O.f (dotted line). It is clear in the figure that a better fit is obtained 
for the template of 2 Myr. Figure 2b shows the plots obtained after fitting a 3 Myr template with E(B — V) = 0.6 (dotted 
line), 4 Myr with E(B ~ V) = 0.5 (solid line), and 6 Myr with E(B ~ V) = 0.4 (dashed line) to the data of C2. In this case 
it is clear that the best fit is obtained for the 4 Myr template. Figure 2c shows the results of fitting templates of 2 Myr with 
E(B - V) = 0.5 (solid line) and 8 Myr with E(B - V) = 0.1 (dotted line) to the C4 data. The best fit this time is found for 
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PA 160 




Figure 3. Ha continuum image showing the locations of the slits and the extraction apertures. 



from 0.0 to 1.6 and the parametrized galactic extinction law of lSeatonl (I1979I ). Each extracted 
spectrum was m o delled using a combination of a young and an old stellar population (see 



Tadhunter et al. 



(120051 )): we used a 12.5Gyr model for the old stellar population (OSP), 
adding YSP ranging in age from O.OOlGyr to 6Gyr (we will refer to this component as 
young component). For each spectroscopic aperture the continuum flux was measured in 
several wavelength bins (~50) chosen to be as evenly distributed in wavelength as possible, 
and to avoid strong emission lines and atmospheric absorption features. We then fitted these 
data using a minimum x 2 technique to determine the best fitting age and reddening, as well 
as the relative contributions of each of the two populations. The result of the modelling the 
spectrum extracted for aperture Ap4 is shown in Figure 4 as an example. 

Although the spectra have the advantage of considerably better wavelength sampling 
than the HST images, they have the disadvantage that several stellar populations of different 
age/reddening may be included in the large spectroscopic aperture. This can lead to an age- 
reddening degeneracy in the sense that adequate fits can sometimes be obtained for either a 
small proportion of a relatively young, highly reddened YSP plus an old stellar population, 
or a much larger proportion of an intermediate age, but less reddened YSP. Figure 5 presents 
the details of the mininmum x 2 e d modelling for aperture Apl and provides a good illustration 
this degeneracy problem (see also Tadhunter et al. 2005), since the two loci that represent 
adequate fits (x 2 e d~ 1) are clearly separated in the figure: at around ~0.1 Gyr with E(B-V) 
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Figure 4. Upper panel: modelling result for the spectum extracted from Ap4 of Fig 3. The green lines represent the young and 
old templates, while the red line represent their sum. The purple boxes represent the bins used for the fit, and the histogram 
on the x-axis represents the modulus of the error measured in each bin. In this case the young component has an age of 0.2 Gyr 
with a reddening of E(B — V) = 0.1. The contribution of this component to the total fiux(red line) is 34.6 % in the normalising 
bin (4600-4700A)and the 

Xj5 e d in this case is 0.4. Lower panel: detailed fit in the wavelength range 3700-4300. The fluxes are 
presented in wavelength units 



= 0.2 and ~1.5 Gyr with E(B-V) = 0.0. To distinguish between groups of templates that 
provide good fits, the model fits to the data were examined in detail, selecting those models 
which provide the best fit to the important absortion features of the spectra (for example 
Call K, G-band, higher order Balmer lines). Figure 6 shows the wavelength range 3700A - 
4300A for two of the best solutions obtained from the SED modelling of Ap3, PA160. It is 
clear from the plot, just by examining the features of the fit for the Call H&K absortion lines, 
that in this case the model comprising 12.5 Gyr OSP plus the 0.04 Gyr YSP component age 
fits the data better than the 1.4 Gyr YSP model. However, there are cases (Ap2, Ap3, ApB) 
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Figure 5. Minimum chi-squared fitting results for Apl in Figure 3. Each contour represents a value for X^e^- Good results for 
the modelling are considered X^ed °^ one or ' ess ( see discussion in Tadhunter et al. 2005). The two loci where good results are 
found arc placed around 0.1 Gyr with E(B — V) = 0.2 and 1.5 Gyr with no reddening 

where even younger (~0.005 Gyr) age YSP — corresponding to the blue A/B supergiant 
phase — provide an adequate fit to the SED. For these latter cases we cannot distinguish 
between YSP with ages ~0.005 Gyr and ~0.04 Gyr using the detailed fits. 

Table 4 sumarizes the results of the modelling for all of the extracted spectroscopic 
apertures. We find ages for the YSP ranging from 0.005 Gyr to 1.5 Gyr with reddenings 
in the range of E(B-V)=0.0 to 1.2. We also measured the mass contribution of the YSPs 
to the overall stellar mass in the aperture, finding results ranging from a negligible mass 
contribution (< 1%) in some apertures, to the YSP contributing a sizeable fraction of the 
total mass in others. The wide range of ages found for some of the apertures reflects the fact 
that good fits can also be achieved for two combinations of age/reddening (the degeneracy 
problem discussed above), and also depends on the nebular continuum subtraction (for Apl 
and Ap2). The reader may notice that the results for NUC AP in Figure 3 are not presented 
in Table 4. After the nebular subtraction, the flux contribution of the YSP in the normalising 
bin is small (< 10%). Therefore we find good fits for the entire range of combinations OSP 
and YSP templates. 

It is notable that for the region to the south of the nucleus that encompasses the SSC 
detected in the HST images (i.e. Ap2, Ap3 and ApB ), the minimum age that provides 
an adequate fit to the SED of the diffuse light (~0.005 Gyr) is consistent with the ages 
of the SSC measured in our photometric study, although older YSP ages (>0.04 Gyr) are 
also acceptable. On the other hand, we can only obtain adequate fits to the SEDs of the 
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Age of 


E(B-V) a 


%YSP 




YSP 




of total 




(Gyr) 




mass 


PA 160 Apl 


0.1 - 0.3 


0.0 - 0.4 


1.0 - 2.0 


Apl-neb 6 


0.2 - 1.2 


0.0 - 0.4 


1 - 69 


Ap2 


0.001 - 0.2 


0.0 - 1.2 


0.05 - 3.5 


Ap2-neb'> 


0.2 - 1.0 


0.3 - 0.7 


1 - 60 


Ap3 


0.005 or 0.04 - 0.2 


0.4 -0.8 or 0.0 - 0.5 


0.1 - 1.5 or 1.5 - 15.0 


Ap4 


0.05 - 0.5 


0.0 - 0.6 


2 - 10 



PA 230 ApA 0.05 -1.5 0.0 - 0.4 2.0 - 90.0 

ApB 0.005 or 0.04 - 1.5 0.0 - 1.0 or 0.0 - 0.7 0.01- 1.2 or 7.0 - 90.0 



Table 4. Results from the analysis of the spectra extracted from the apertures in Figue 3. 

a The large range in reddening should be understood as it follows. For the younger solutions there is a wide range of reddenings 
which gives good results. However, when moving towards older ages the amount of reddening needed for the template to match 
the data decreases as does the range. 

b Two different sets of results are presented for the apertures 1 and 2. One is obtained from modelling without subtracting any 
nebular continuum from the spectra, while the other (labelled Apl-neb and Ap2-neb) is obtained after having subtracted the 
maximum nebular continuum. Thus the true results are likely to fall between the the resuls presented in the table 

other extended apertures (Apl, ApA Sz Ap4) using older YSP with ages in the range 0.04 
- 1.0 Gyr, and we cannot rule out the idea that YSP sampled by the diffuse light have the 
same age (~0.1 - 0.3 Gyr) for all apertures. 

Overall, the range of YSP ages we have determined for the diffuse light from our spectra 
of PKS1345+12 overlaps with the 1 - 300 Myr range determined for the YSP in the spec- 
troscopic study of a sample of transitional ULIRG QSOs by Canalizo et al. (2001); it is also 
consistent with the 0.5 - 1.5 Gyr age determined by Tadhunter et al. (2005) for a single 
extended aperture to the SE of the nucleus of PKS1345+12 itself, although this latter study 
considered only unreddened YSP. 

Finally we emphasise that the diffuse light contains the dominant YSP component in 
terms of mass. A lower limit on the total mass in the diffuse YSP component integrated 
across all our spectroscopic apertures is Mj^J > 1.4 x 10 9 M Q , and may be considerably 
larger if the regions of the galaxy not sampled by the slits are included. In contrast, the total 
mass in the three SSC with accurate age estimates is orders of magnitude less (8 x 10 6 < 
Ml 3 / < 2 x 10 7 M ). 



3.3 Gas Kinematics in the Halo 

The emission line kinematics of PKS1345+12 have been described in detail by Holt et 
al. (2003). As well as a quiescent emission line halo with relativ ely small r a dial v elocity 



amplitudes (<400 km s" 1 ) and linewidths (<500 km s" 1 FWHM UHolt et all ((20031) found 
extreme kinematic components at the location of the western nucleus with velocity shifts up 
to ~2000 km s _1 and line widths of ~2000 km s _1 (FWHM) providing evidence for warm 



16 J.Rodriguez Zaurin, J.Holt, C.N.Tadhunter and R.M. Gonzalez Delg ado 

Balmer lines 

1.5«10 



l.O'lO 17 H 



5.0«10 1S - 




3700 



3800 



3900 4000 4100 

Rest wavelength 



4200 



4300 



1.5«10 



1.0* 10 17 r 



5.0-10 18 |- 



3700 



Balmer lines 




3800 



3900 4000 

Rest wavelength 



4100 



4200 



4300 



Figure 6. Detailed fits for the spectrum extracted from Ap3. Upper panel: detailed fit for a model with a young stellar 
component of 0.04 Gyr (E(B-V)=0.2). Lower panel: detailed fit for a model with a "young" stellar component of 1.4 Gyr (no 
reddening). The fluxes arc presented in wavelength units. It is clear from the figure that the model with the younger YSP 
(upper panel) fits the absortion features best, despite the fact that both models provide adequate fits to the overall SEDs. 



gas outflows in the narrow line region. As mentioned in Section 3.1, we detect emission lines 
at the locations of the young star clusters visible in our HST in both the images and the 
long slit spectra. Therefore it is interesting to compare the kinematics of the star clusters 
with those of the other kinematic components in the galaxy. 
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Figure 7. Spatial variations of the [O n] and Ha radial velocities of the extended gaseous halo along PA 160. [O n] (open 
diamonds), Ho (open diamonds + open squares). The err or bars are sm all and are incorporated in the points. Overplotted is 
the radial velocity of the Hi 21cm absorption (filled circle; Mirabel (1989)), the three components of [Oin] in the nucleus (open 
triangles) and the extended emission line regions: the extended [Om] emission in the region to the NW (large open diamond) 
and the extended Ha emission in the line splitting region to the SE coinciding with the position of the cluster CI (large filled 
diamonds). Note: Positive kpc = south, negative kpc = north. The vertical dotted line represents the centroid of our nuclear 
aperture. The horizontal dotted line represents the velocity at the position of the centroid. 

Figures 7 and 8 compare the emission line kinematics at the locations of the clusters 
with those of the spatially resolved diffuse gas detected in the Ha and [OII]A3727 emission 
lines in the halo of the galaxy along PA 160 and PA 230. Intriguingly, we find clear evidence 
for line splitting in both Ha and [Nil] at the locations of clusters C1/C2 and C4. The line 
splitting is clearly visible in the extracted spectra of these regions shown in Figures 9 and 10 
and also the grey scale representation of the long-slit spectra shown in Figure 11. Important 
features of the line splitting regions include the following. 

• Linewidths. In both regions there is a relatively broad kinematic component (320 < 
FWHM < 400 km s _1 ) along with a narrower kinematic component (FWHM < 150 km 
s _1 ) that is spectrally unresolved. 

• Radial velocity shifts. In both regions the broader component has a radial velocity 
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and line width consistent with that of the diffuse gas at larger and smaller radial distances 
from the nucleus, whereas the narrower component is significantly shifted relative to the 
diffuse gas: blueshifted by ~450 km s _1 for region C1/C2, and redshiftedby ~300 km s _1 for 
region C4. 

• Ionization state. Whereas the narrow kinematic component has relatively small [NII]6584/Ho; 
and [OIII]/H/3 ratios ([NII]6584/Ha < 0.5 and [OIII]/H/3 < 1), consistent with photoioniza- 
tion by hot stars in an HII region (see Veilleux et al. 1999), the broader component has line 
ratios ([NII]6584/Ha > 0.7 and [OIIIj/H/3 < 3) that are more consistent with LINER nuclei 
and the diffuse gas at other locations in the host galaxy. Overall the LINER-like line ratios 
for the broader diffuse gas component indicate ionization by shocks, or photoionization by 
an AGN at low ionization parameter. 



All of these features are consistent with the idea that the narrow emission line components 
represent HII regions associated with the young star clusters, and that these SSC are moving 
at high velocity relative to both the ambient gas and the rest frame of the galaxy defined 
by the narrow emission lines and the HI 21cm absorption line detected in the nucleus. In 
this case, we are observing the clusters before the phase in which th e s upernovae associated 



with the massive stars in the clusters eject the ISM ((Goodwin, 



1997allbl: Bastian fc Goodwi: 



20061 ). Indeed, given that the supernova phase is expected to start ~3 - 4 Myr after the birth 



of a star cluster (Leitherer et al. 1999), the ages we have determined for the SSC ( <6 Myr) 
are consistent with such a supposition. Given their measured radial velocities and estimated 
ages, the clusters must have formed within a few kpc of their current locations, unless their 
tangential velocities are considerably larger than their radial velocities. This latter feature 
rules out the idea that the clusters formed at much larger radii in the outer halo of the 
merging system, then fell into the nuclear regions; the gas out of which the SSC formed must 
have been moving rapidly prior to the formation of the star clusters. 

The alternative explanation, that we have caught the SSC in the ejection phase, and that 
the line sp l itting is caused by the ejection of ga s by SSC formed in the quiescent, diffuse gas 



( iGoodwinl . Il997al lbl: iBastian fc Goodwin! 



2006) is difficult to reconcile with the observations. 



In particular, in such a case it is not clear why the component apparently at rest relative to 
the ambient gas has relatively broad lines and a [NII]6548/Ha ratio consistent with shocks or 
AGN photoionization, whereas the component shifted relative to the rest frame (presumably 
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Figure 8. Same figure than 8 but for PA 230. [O II] (asterisks), Ha (asterisks-l- open squares). The err or bars a r e sma ll and 
are incorporated in the points. Overplotted is the radial velocity of the Hi 21cm absortion (filled circle; Mirabel (1989)), the 
three components of [O III] in the nucleus (open triangles) and the extended emission line regions: the extended [O III] emission 
(large asterisk) and the extended Ha emission in the line splitting region to the SW coinciding with the position of the cluster 
C4 (large stars). Note: Positive kpc = south-west, negative kpc = north-east. The vertical dotted line represents the centroid 
of our nuclear aperture. The horizontal dotted line represents the velocity at the position of the centroid. 



the ejected component in this scenario) has narrow lines and an HII region-like [Nil] /Ha 
ratio characteristic of stellar photoionization; one would expect the reverse to be the case. 

Finally we note that, despite the apparently "free-floating" status indicated by their 
emission line kinematics, the SSC/HII regions are linked with a more extensive interstellar 
medium (ISM) in the host galaxies. First, there is clear morphological evidence from our 
high resolution ACS images that the young star clusters are associated with dust lanes in the 
host galaxy (see Figure 1). Second, at the location of C1/C2, the greyscale representation 
of the long-slit spectrum (see Figure 11) shows a clear enhancement in the flux of both the 
broad and narrow components, suggesting a link between the HII regions and the more 
diffuse gas. 
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Figure 9. Ha model of the region exhibitting line splitting coinciding with the position of CI. The faint line traces the observed 
spectra and the bold line represents the overall model profile. The six components, two for each line, are also overplotted — 
the dotted line traces the broader component consistent with the kinematics of the emission line gas at other locations in the 
galaxy, and the dot-dashed line traces the narrow component blueshifted by ~450 km s _1 with respect to the ambient gas. 



4 DISCUSSION 



4.1 Age determination for the YSPs in ULIRGs 

It is inte resting to compare our results for PKS1345+12 with previous studies of YSP in 



ULIRGs. 



Surace et al 



(119981 ) carried out an HST imaging study concentrating on bright 



star forming knots for a sample of nine "warm" ULIRGs, including PKS1345+12, using 



high resolution B- and I-band images taken with the Wide Fie 



HST. The stellar synthesis models of 



Bruzual fc Chariot 



d Planetary Camera on the 



(119931 ) were used to estimate the 



ages and masses of the bright knots in colour-magnitude diagrams. However, It is difficult to 
distinguish between stellar age and redd ening using only B- and I-band observations. Thus 



the ages presented in 



Surace et al. 



(119981 ) are upper limits and, as discussed in that paper, the 
effect of correcting for the reddening could substantially modify the age and mass estimates, 
resulting in younger ages and smaller masses. For the whole sample, they derive a median 
upper limit for the age of ~ 3 x 10 8 yr and masses in the range 10 5 < M c \ usteTS < 10 9 M . 
In the case of PKS1345+12 itself, they found ages ~ 10 7 < t ssc < 10 8 yr and masses 



~ 10 7 < M, 



clusters 



< 10 8 M . Since we have measured significant amounts of reddening 
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Figure 10. Ha model of the region exhibitting line splitting coinciding with the position of C4. The faint line traces the observed 
spectra and the bold line represents the overall model profile. The six components, two for each line, are also overplotted — 
the dotted line traces the broader component, and the dot-dashed line trace the second narrow component redshifted by ~300 
km s — 1 with respect to the spatially extended narrow component 




Figure 11. 2-D optical spectra for PA 160(left) and PA 230(right), centred on the redshifted Ha+[NII] emission lines. The x 
axis is the wavelength direction increasing to the right, and the vertical axis is the spatial direction with SE to the top and 
NW to the bottom. The images show the regions exhibiting Ha line spitting, coincident with the SSC, and demonstrate the 
flux enhancement of both kinematic in the line splitting region in PA 160 (at the location of C1/C2). 



for the knots in PKS134 5+12, this ga l axy is clearly a case in which the reddening affects 



the results presented in 



Surace et al. 



(119981 ) . As expected, after taking into account the 



reddening, we deduce younger ages and lower masses for the SSC in PKS1345+12. 

It is well known that large amounts of gas and dust surrou n d the AGN in the western 



nucleus of 



Holt et al. 



KS1345+12 (Evans et al. 



1999 



Surace fc Sanders! . Il999l : IScoville et al. 



2000; 



20031 ). Our new results demonstrate the importance of dust on a larger, ~5 - 



15 kpc scale, in the halo of the galaxy. Not only do the SED fits require a reddening in the 
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range 0.2 < E(B — V) < 0.5, but there are dust structures in our images coincident with 
the SSC. Therefore it is dangerous to use single colour-magnitude diagrams to determine 
the ages of the YSPs in such objects. 

As well as the importance of taking into account reddening, our results also demonstrate 
that imaging and spectroscopic data provide complementary information about the YSPs; 
clearly a combination of the two types of data is required to make a full census of the YSPs 
present in systems of this type. Whereas on the basis of photometric analysis we find robust 
ages <6 Myr for the three SSC with good SEDs, our spectroscopic analysis places a lower 
limit of 40Myr on the ages of the YSP associated with the diffuse light in three of the 
regions sampled by our spectra. Therefore, our results demonstrate that the highest surface 
brightness regions (in this case the SSC) can give a misleading impression of the properties 
of the YSP component associated with the diffuse light in the galaxy as a whole. 



4.2 The past and future of star formation in PKS1345+12 



Combining both imaging and spectroscopic techniques, we find evidence for two distinct 
phases of recent star formation activity in PKS1345+12: one occuring less than 6 Myr ago 
that is associated with the SSCs and perhaps also linked to current the AGN and ULIRG 
activity; the other occuring more than 40 Myr ago and associated with the diffuse light. 
This is consistent with previous HST imaging studies which have shown evidence for more 



than one episode of star formation a ctivity in oth e r merging systems inc 



(Whitmore et al. 



(IWilson et al 



1999|), NGC7252 (Miller et al. 



1997 



Maraston et al. 



uding the Antennae 



20011 ). and Arp220 



20061 ). Therefore our results fit in with the emerging trend that the star 



formation in merging systems is complex and multi-modal. 

At this stage it is important to add the caveat that, although we have assumed that all 
of the YSP in PKS1345+12 formed in starbursts triggered by the merger, an alternative 
explanation is that the older YSP component associated with the diffuse light represents 
the captured disk of one of the merging galaxies i.e. the diffuse light YSP is not directly 
associated with the merging process. To test this latter scenario we have attempted to model 
the spectra of the diffuse light for all ape rtures using templat e spectra for disk galaxies of Sa 



and Sb morphological types taken from 



Kinney et al. 



(11991 . We find good fits with either 



Sa or Sb templates for Ap 1, Ap2, a nd Ap A and Ap B. Although no good fits are found 



(Xred > 2) for the other apertures, the 



Kinney et al. 



(119961 ) templates represent averages for 
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several galaxies of the same morphological type, and each galaxy type template encompasses 
a range of individual galaxy spectra, some of which deviate substantially from the average 
SED. Therefore we cannot absolutely rule out the idea that the YSP detected in the diffuse 
light for these apertures are associated with the disrupted disk(s) of one or more of the 
merging galaxies. Overall it is entirely plausible that the YSP associated with the diffuse 
light represent the captured stellar population from one or more merging disk galaxies, rather 
than stars formed in the merger itself 

Are our results consistent with the numerical simulations of major galaxy mergers? For 
the case of a merger between two galaxies with substantial bulges — the likely scenario for 
PKS1345+12 — the simulations predict that the major merger-induced starburst occurs as 
the two nuclei finally merge together, but that there will also be a lower level of interaction- 
induced star formation activity at an earlier stage due to the tidal effects of the interaction on 



the d isks of the individual merging galaxies (IBarnes fc Hernquist 



1996 



Mihos &: Hernquistl . 



19961 ) . Therefore it is plausible that the younger YSP represented by the SSC and the ULIRG 



activity are all related to the high gas densities produced in the final stages of the merger 
as the nuclei merge together. On the other hand, assuming that they are not captured disk 
populations, the older YSP associated with the diffuse light may represent the earlier phase 
of lower level interaction-induced star formation. In this case, the models predict that we 
should expect an even higher level of star formation a nd AGN activity over the ~ 10 7 yr it 
will take for the nuclei to coalesce (jEvans et al.l (119991 )). 



4.3 Order of events in triggering the AGN 



Several recent studies of YSP have been used to estimate the timescales and order of 
eve nts for triggering the ac tivity in radio galaxies. For their sample of three radio galax- 



ies. 



Tadhunter et al. 



20051 ) found relatively old post-starburst ag es (~0.3 - 2.5Gyr) com 



pared with the lifetimes of the radio sources. 



Emonts et al. 



(120061 ) also found a significant 



(~0.3 Gyr) delay between the starburst event and Radio-AGN activity for the radio galaxy 
B2 0648+27. These results are consistent with a scenario in which the radio activity is trig- 
gered rel atively late in the merge r sequence, following the major merger-induced starburst. 



However, 



Tadhunter et al. 



(120051 ) also discussed the possibility tha t some radio sources have 



undergone multiple epochs of activity. Indeed, 



Stanghellini et al 



(120051 ) detected low sur- 



face brightness extended radio emission surrounding several objects in their sample of oth- 



24 J.Rodriguez Zaurin, J.Holt, C.N.Tadhunter and R.M. Gonzalez Delgado 
erwise compact {D < 1 kpc) radio sources, including PKS1345+12. Moreover, this galaxy 
shows some evidence for continuity between the compact and extended emission, which, 
if confirmed, would suggest that PKS1345+12 has had AGN- Radio activity over a longer 
timescale than the estimated age of the GPS source (<0.1Myr). Given the evidence for 
both multiple epochs of star formation, including some YSP with ages less than 6 Myr, and 
multiple epochs of radio jet activity, we conclude that there is no clear evidence for a delay 
between the merger-induced starburst and the start of the AGN phase in PKS1345+12. 



4.4 SSC kinematics 



In Section 3.3 we presented evidence that the SSC and their associated HII regions are 
moving at high radial velocities relative to the diffuse ISM detected (in projection) at similar 
locations in the halo of the host galaxy, despite the morphological evidence for associations 
with diffuse ISM components. It is interesting to consider how such rapidly moving star 
clusters might form. 

The first possibility recognises the fact that major galaxy mergers are violent events, 
and considerable dissipation must take place before the gas finally settles into a stable con- 
figuration in the merged system. At all stages in a merger, but particularly the final stages 
before the nuclei of the two merging galaxies coalesce, it is likely that rapidly moving gas 
streams/tidal tails will rain down on the central regions of the host galaxy. As individual gas 
streams collide with the am bient gas or o ther gas streams, they will generate shocks that 
may trigger star formation (JBarnesJ, 12004 ). The ionizing photons generated in the shocks 
as the tidal streams interact with the ambient gas may also lead to a local enhancement 
in the emission line luminosity of the ambient gas, perhaps explaining the effect seen in 
Figure 11. Therefore the formation of the SSC detected in PKS1345+12 may be associ- 
ated with the overall settling/dissipation process that must accompany the evolution of the 
ISM following a major galaxy merger; it is possible that we are witnessing the events that 
will eventually lead t o the f ormation of a high velo city dispersion globular cluster systems 



(e.g 



Ashman fc Zepj (119921 ) : IZepf &: Ashmanl (119931 )) as the host galaxy settles down to its 



final state quiescent elliptical galaxy. 

A second possibility is suggested by the detection of diffuse radio emission extending ~35 
arcseconds (~83 kpc) to the north, and ~25 arecseconds (~60 kpc) to the south, of the nu- 
cleus of PKS1345+12, and encompassing the region of the SSC (jStanghellini et all 120051 ). As 
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discussed in the previous section, although the major, high surfac e brightness radio activit y 



is now taking place on a sub-kpc scale close to the western nucleus ( iStanghellini et al. 



200,4 . 



the presence of diffuse radio emission on a much larger scale provides evidence that there 
was an earlier phase of activity in which the relativistic plasma jets escaped into the halo of 
the merging system. In this phase, the radially expanding jets and laterally expanding jet 



cocoon would have interacted strongly with the warm/cool 



shock ing the gas and possibly triggering star formation (IRees 



SM associa t ed with the merger 



1989 



Begelman fe Cioffil . 



19891 ). Because the shocked gas is accelerated to high velocities by the jet-cloud interactions, 



it would be expected that any SSC formed in such interactions would be moving at a high 
velocity relative to the ambient ISM, as observed in PKS1345+12. Moreover, interaction 
with the expanding cocoon could produce both positive and negative radial velocities, de- 
pending on the interaction geometry. Therefore, the fact that C1/C2 is blueshifted, while 
C4 is redshifted, is not a problem for this model. 

Note that the importance of jet-induced star formation as a mechanism remains con- 
troversial. The direct observational evidence for the mec hanism remains sparse, confined to 



a few objec t s obse rved across a large range of redshifts (IBest et al. 



1997 



Croft et al 



Dey et al. 



1997 ; 



20061 ). There is also a recognition that hydrodynamic effects of the jet-cloud 



interactions, as well as the heating effects of the accompanying quasar activity, can be as 
destructive as they are conducive to star formation. However, the recent d etection of high 



veloc it y neutral outflows ag ainst the cores of several compact radio sources ( IMorganti et al 



2005l a. iMorganti et al.ll2005l b) — including PKS1345+12 itself — demonstrates that, despite 



being accelerated to high velocities, the gas behind a shock can at least cool to a neutral 
phase. This lends weight to the plausibility of jet-induced star formation as a mechanism. 



5 CONCLUSION 

We have reported a detailed investigation of the YSP in the ULIRG PKS1345+12, combining 
information from both spectroscopic and imaging methods. Our key conclusions can be 
summarised as follows: 

• Reddening. We find that the three SSC with good photometric information in our 
study are all significantly reddened. Clearly it is important to have sufficient photomet- 
ric information with wide spectral coverage to remove any age-reddening degeneracy and 
determine accurate YSP properties, even for clusters in the extended halo of a galaxy. 
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• High resolution imaging versus spectroscopy. This study provides a clear demon- 
stration of the fact that high resolution HST imaging studies tend to be biased towards 
the youngest and/or most massive individual star forming clusters, and may give a mislead- 
ing impression of the star formation histories of the dominant YSP in terms of total mass 
(sampled by the diffuse light in PKS1345+12). 

• The star formation history of PKS1345+12. We find evidence for a complex 
star formation history of the YSP in PKS1345+12, with at least two major phases of star 
formation, one of which may be associated with star formation in the merging disk galaxies 
at an earlier stage of the merger. 

• Order of events. Unlike some other radio galaxies, we find no clear evidence in 
PKS1345+12 for a time delay between the major merger-induced starburst and the AGN 
activity. 

• Star cluster kinematics. On the basis of our analysis of the emission line kinematics, 
we deduce that the SSC are moving at high velocity relative to ambient gas and the rest 
frame of the host galaxy. The extreme SSC kinematics may be related to the process of star 
formation as infalling gas settles in the nuclear regions in the final stages of the merger, or 
to star formation triggered by the expanding radio source as it drives shocks into the ISM. 
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